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Abstract 

We propose a new method for measuring CP violation in neutrino oscillation 
experiments. The idea is to isolate the term due to the CP-violating phase 
out of the oscillation probability by taking difference between yields of two (or 
three) detectors at path-lengths L = 250 ssCev ) (lo^Tv^) 

L/3 (and also at 2L/3 in the case of three detectors). We use possible hierar¬ 
chies in neutrino masses suggested by the astrophysical and the cosmological 
observations to motivate the idea and to examine how the method works. 
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CP violation in the lepton sector is an unexplored fascinating subject in particle physics. 
If observed, it should shed light on the deep relationship between quarks and leptons, the 
most fundamental structure of matter which we know to date. Moreover, it is suggested 
that CP violation in the lepton sector is the one of the key ingredients of the mechanism 
for generating baryon number asymmetry in the universe []^. 

A viable way of observing CP violation in the lepton sector is to utilize the phenomenon 
of neutrino oscillation. It was pointed out in refs. that the difference between oscillation 
probabilities of the neutrinos and its antineutrinos is proportional to the leptonic analogue of 
the Jarlskog factor |^, the unique (in three-flavor neutrinos) phase-convention independent 
measure for CP violation. Recently, measuring CP violation in long-baseline neutrino 
oscillation experiments is of considerable interest in the literature 

However, there exists a potential obstacle in measuring CP violation in long-baseline 
neutrino experiments. It is the problem of the contamination due to the matter effect. 
Since the earth matter is not CP symmetric its effect inevitably produces the fake CP 
violation []^ which acts as contamination to the genuine CP-violating effect due to the 
leptonic Kobayashi-Maskawa phase |n|. Even worse, the matter effect dominates over the 
CP phase effect in a certain region of the mixing parameters in the —>■ z/g experiment 


In this paper, we suggest a novel way of measuring CP violation in neutrino oscillation 
experiments by proposing the multiple detector difference method.[] We will show that our 
method is relatively free from the problem of the matter effect contamination in long-baseline 


^Preliminary descriptions of the two-detector difference method were given in ref. |T^. We should 
mention that the idea of placing multiple detectors in long-baseline experiments is not new. For 
example, it appeared in Brookhaven proposal [^. But their motivation and the basic idea be¬ 
hind the use of the multiple detectors is entirely different from ours, and they do not discuss the 
possibility of measuring CP violation. 
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neutrino oscillation experiments. 

One way of avoiding the problem of “matter effect pollution” is to look for the oscillation 
channel in which genuine CP violating effect dominates over the matter effect. This idea is 
examined in detail in ref. p| within the restriction of the neutrino mass hierarchy to that 
motivated by hot dark matter and the atmospheric neutrino anomaly. It is found that under 
the constraints from the terrestrial experiments the unique case where the CP phase effect 
dominates is the Ve channel in the region of large-sis and arbitrary-S23 (the region (B) 
to be defined later). The Vr channel is relatively free from matter effect contaminations 
but they are not negligible. Unfortunately, the expected CP violating effect is at most ~ 
1% in this type of mass hierarchy due to the strong constraints on mixing angles from the 
terrestrial experiments. 

These discussions of the absolute and the relative magnitudes oi CP violation is based 
on the V — V difference method.^ The major experimental problem with this method is the 
difficulty in determining relative normalization of the neutrino and the antineutrino beams. 
If the CP violation to be measured is of the order of a few % it would require to calibrate the 
flux of neutrino beams to the accuracy better than it, which would be extremely difficult, if 
not impossible, experimentally. 

The multiple detector difference method which we shall discuss in this paper aims to 
overcome this problem. Since the absolute flux of the neutrino beam is hard to determine 


to the accuracy better than ~ 10% |T^ we may have to give up the comparison between 
the two beams, or two different experiments, if we want to do measurement with a few % 
level accuracy. Therefore, we stick to use a single neutrino beam, for example, in the 
experiment. Then, how can one measure CP violation to such precision, or even to the 
accuracy of 1% level? 

In this paper we confine ourselves into the three-flavor mixing scheme of neutrinos. To 


^Related but different proposals in this context have been discussed in ref. Q. 
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develop the multiple detector difference method we work with the mass hierarchy 


AM^ = Amfg ~ Am23 S> Am^2 = Am^, 


( 1 ) 


where Am^^ = {i,j = (1,3), (2,3), (1,2)), motivated by the solar and the atmospheric 

neutrino observation [^,|T^ and the neutrinos as hot dark matter in the mixed dark matter 
cosmology . 


We hrst focus on the case of mass hierarchy motivated by the dark-matter and the 
atmospheric neutrino observation, AM^ = 5 — 100 eV^ and Am^ = 10“^ — 10“^eV^. To hnd 
a hint on how to isolate the CP violating term we hrst ignore the matter effect and analyze 
the structure of the neutrino oscillation probability in vacuum. With the mass hierarchy the 
oscillation probability in vacuum for the long-baseline experiments can be written as 


P{i >0 z/«) = Alia + 5/3a(l - COS A) C^a sin A 


( 2 ) 


with 


A = 


/S.m?L 
2E ’ 


where L denotes the path-length of the baseline, and B^a and Cjia are constants which 
depend on the mixing angles and the CP phase. We note that Cpa = 2J, up to the sign, 
where J indicates the Jarlskog factor whose explicit expression will be given later. (See 
eq.(|T^.) The rapid oscillation due to large AM^, 

L \ f E \-^ 


AM‘^L „/AM2\ 
= 127 


AE 


/ AMM / L \/E\-^ 
l^leV^ ) UOOkmj UGevj ’ 


( 3 ) 


is averaged out which produces the hrst term in (|^). 

We want to isolate the last term, J-term, which is the measure for CP violation from 
the others. It is a simple matter to observe that the best way to carry this out is to do the 


measurements at A = | and |7r: 


dvr. 


^ ^a-j A ^ ) Alia T BfiQ 

Tl 

P{yi3 ^ A —) Apa T B^q 


2J 

2J 


( 4 ) 
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where we took a particular sign for the J-term. Therefore, the difference AP between the 
oscillation probabilities at A = |7r and at A = | is nothing but the CP violation 4J. 

Of course, A is a function of L/E. Therefore, the measurement of the difference P(A = 
— P(A = ^) = AP 2 can be done either by varying L or P, or both.^ However, we 
argue that a measurement of OP violation with accuracy better than a few % would require 
measurements of AP 2 by placing two detectors at A = | and A = |7r. If we try to measure 
AP 2 by varying the energy of the neutrino beam with the single detector the energy, of 
course, have to be varied by factor of 3. By such re-adjustment of the neutrino beam energy 
the relative normalization of the beam would become uncertain to the order of ~10%. 
Therefore, the best thinkable way of avoiding the uncertainty of relative normalization of 
the neutrino flux is to do measurement at 2 detectors, using the same neutrino beam, one at 
A = Itt, and the other at A = |. For the KEK-PS—s>Superkamiokande experiment in which 
L =250 km, the neutrino beam energy should be tuned to E = 1.35( ^Q^^y2 )~^ GeV so that 
the location of Superkamiokande just corresponds to A = |7r. For the MINOS experiment 
with L = 730 km the beam energy to be used is P = 3.94 ( ^q^^^^ 2 )~^ GeV. Then, the 
second detector to be build should be located at 1/3 of the baseline, L 2 = 83.3 km for the 
KEK-PS—s>Superkamiokande and L 2 = 243 km for the MINOS experiments.^ 

Of course, one cannot make neutrino beam so monochromatic; it must have spread in 
energy. However, it seems possible to make the neutrino beam spread as small as ~ 20 % 


of the beam energy [Q. Therefore, at least it is worth to think about the possibility. 

Does the matter effect give rise to the serious contamination to the genuine CP violating 
effect in the 2 detector difference method? It appears that the problem has a better shape 


^The possibility of observing CP violation in neutrino oscillation experiments by measurement at 
differing path-length, or by varying energy has been pointed out in ref. [|^]. 

^Their geographical locations are, respectively, at around the city of Honjo in Saitama prefecture 
and a midpoint between Westfield and Packwaukee, about 55 Miles north of Madison, Wisconsin. 


5 








compared to that in the u — U difference method. To understand this point we write down 
the expression of neutrino oscillation probability with the correction of the earth matter 
effect. The expression is based on the adiabatic approximation and is valid to first order 
in matter perturbation theory. If we define the neutrino mixing matrix as z/q = Uai^i the 
oscillation probability is given by 




= 


-2 ^ Re[U^,U:,U;M + MUUU6V)^p,,s 
1=1,2 f 

1 




-2J 


f Am? 


sm 


sm 


AE 

2 


L]+^aL[\Ue2\' -\U, 


+ (^L[\Ue2V - \UeiVj cos(^ 


2 

/ Am^ 


>el\ 


sm 


f Am? 


-4Re(t/f/f/5I/)„/3;i2sin 


, / Am^ ^ 
L 


2E 


-2lm{UUU6V)^p,i2 sin(^ 


AE 
f Am? 


2E 


(5) 


where {UUU6V)a^ -ij represent first-order corrections due to the matter effect and their 
expressions are given in ref. |§|. More precisely speaking, we made the following approxi¬ 
mation to derive eq. (|]): We took average over rapid oscillations with the period (j^), which 
produces the first two terms in (||). We ignored the terms of the order of -^2 because of 
the extreme hierarchy between the dark matter mass scale and the matter potential, 


Ea 

AM2 


= 1.04 X 10 


-4 


P 


f E \ /AM^Y 

UGevj I leV ') 


( 6 ) 


^2.72gcm“3 

We used the constant matter density approximation and ignored the terms of order (aL)^ 
or higher, where L is the path length of the baseline and a = \/2GFNe with being the 
electron number density. We note that 


P 


(- 


L \ 


(7) 


aL = 0.132 -^-- I (-^— I 

V2.72gcm-3 ) V250kmy 

Therefore, ignoring (aL)^ term should give a good approximation at least for the 

KEK—>Superkamiokande experiment. 

The difference between probabilities AP 2 to be measured at 2 detectors up to the first 
order in matter potential a is given by 


6 





















( 8 ) 


3 TT 

AF 2 (z ^/3 ^ A = -tt) - P{vf} ^ z/„; A = -) 

= 4J - M^^Re[UaiU*^2UliUp2] cos26'i2Ci3 

Ea „ n 
+8-/^ COs2012Ci3 

where (a, /3) = (e, /i), (/i, r) and (r, e). We have used the standard form of the CKM matrix 




C12C13 

-512O13 

5130 *'^ 


— -512023 

- 012 - 523 - 5136 *'^ 

O12O23 - -512523-5130*'^ 

523O13 

( 9 ) 

- 512-523 

- 012023-5130*^ 

-O12S23 - 5120235130*'^ 

O23C13 



for the neutrino mixing matrix. With this parametrization, J is given as 


J = lm[UaiU^2i^^lV/32] ^ ±Ci2Sl2C23S23Ci3Sl3Sill5, 


( 10 ) 


where + sign is for cyclic permutations, i.e., (a,/?) = (e,p), (p., r), (r, e) and — is for anti- 
cyclic ones. 

We note that in the mass hierarchy with dark matter mass scale the reactor and the 
accelerator experiments constrain the mixing angles into three regions on the plane spanned 
by S 33 and S 23 [jl^- Since our present analysis is motivated by the atmospheric neutrino 
anomaly we restrict ourselves into the two regions (A) or (B) : 


(A) small-si 3 and small-S 23 


(B) large-si 3 and arbitrary S 23 

In long-baseline experiments, we expect large (~ order 1 ) oscillation probabilities for 
and —> Ut- channels in the regions (A) and (B), respectively. For brevity we 

shall call these channels the dominant channels in the respective parameter regions, and 
the alternative ones, i.e., —>■ Ur in the regions (A) etc. the minor channels. Note that 

— 4:Re[UaiU*2U^iUi32] is nothing but the coefficient of sin^(^^^L) term in the oscillation 
probabilities, and therefore is of order unity for the dominant channels. For AM^ = 5 eV^, 
4J is at most ~ 10“^ in both regions (A) and (B). See Fig. 1 of ref. [^. We note that 
~ 10 “^ in the mass hierarchy with which we are dealing. 
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Now we can roughly estimate the relative magnitude of the matter and the CP violating 
effects in AP 2 . It depends upon the regions (A) and (B); In the region (A), the CP phase 
and the matter effects are, roughly speaking, comparable unless cos 26^12 is tuned to be 
small. We have to do the experiments in the minor channel —>■ Ur to avoid matter effect 
contamination. In the region (B) on the other hand, C 13 is small, ~ 10“^, and therefore 
the CP phase effect is always dominant over matter effect. One can do experiments in the 
dominant channel in the region (B). However, to do subtraction between the yields of the 
end and the intermediate detectors it may be better to do experiments always in the minor 
channel; One has to subtract a large number from a comparable large number to have a 
small number in the case of dominant channels. 

A special attention has been paid to the long-baseline neutrino experiment because with 
the dark matter motivated mass hierarchy the effect of OP violation cannot be observed 
in short-baseline experiments. It can be shown that CP violating effect is suppressed by a 
factor of 

Let us conhrm by numerical computation that the qualitative results we have obtained 
so far are correct. First we pick up the following two sets of parameters (a) and (b) from 
the allowed regions (A) and (B), respectively. 


(a) sl^ = 3.0 X 10"^ = 2.0 x lO’^ 

(b) = 2.0 X 10-2, s^g = 0.98. 


( 11 ) 

( 12 ) 


They are chosen so that the largest CP violating effect is expected in each region (A) or 
(B), and the same set of parameters are used in our analysis done in ref. §. 

In Fig. 1 we present AP 2 (z/^ —>• z/g) and AP 2 (z/^ —>• Pr) for these two parameter sets 
for the KEK—Superkamiokande distance L = 250 km. We have carried out the calculation 
using the exact solutions obtained in ref. for a constant matter density and we used p = 
2.72 gcm”^ with the electron fraction W = 0.5. We took average over the rapid oscillations 
due to the dark matter scale AM^. In the same hgure we also plot the values obtained by 
our analytic formula (H) to indicate that it gives a reasonably good approximation. 
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In Fig. 2 we present the same quantities but for L = 730 km, i.e., Fermilab-Soudan2 


detector distance (MINOS experiment). This is important for > Uj. channel because r 
cannot be produced with such low neutrino energy as E = 1.35 GeV. The nice feature of 
AP 2 in Figs. 1 and 2 is that the matter effect contamination is relatively smaller than that 
in the u — 9 difference method discussed in ref. [Q. This is particularly true for L = 250 km. 
Therefore, if we can measure AP 2 to the accuracy of ~ 1 % level it is, in principle, possible 
to observe CP violation in neutrino oscillation experiments. 

Now we turn to the question of how the multiple detector difference method works for 
the neutrino mass hierarchy motivated by the atmospheric [0 and the solar |T^ neutrino 
observations, AM^ = 10“^ — 10“^eV^ and Am^ = 10“® — 10“‘^eV^. In this case we cannot 
use the matter perturbation theory developed in ref. [§ because -^2 — 1 — 10 ^ cannot be 
used as an expansion parameter. In this paper we rely on the approximate formula derived 
by Arafune, Koike and Sato 0 who use as expansion parameters aL in and 


AtiPL 


= 6.4 X 10 


-2 


Am^ \ / L 


2E l^lO-W^yl l250km) (iGeVy* 

The general structure of the oscillation probability can be expressed, to leading order 
of these expansion parameters, as 


E \ 


-1 


(13) 


P{.^I 3 9a) = A/3q,(1 — COS A) + i?^Q,Asin A + Cj^aAil — cos A) (14) 

where A = ■ The coefficient of the last term is given by The coefficient 

Bf^a contains the terms proportional to either or and they have similar order 
of magnitude as the last term. Enriched by three terms with different A dependence and 
with similar magnitudes it is unlikely that one can separate the last term by using only 
2 detectors. Then, we have to generalize our method to that of 3 detectors, the ones at 
A = 7r/2, vr, and 37r/2. Then, it is simple to show that 

Q TT /\t77 ^ 

AP 3 (^/. - ^a) ^ P{A = -TT) + 3P(A = -) - 2P(A = vr) = 27rJ^ . (15) 

To have a feeling of the magnitude of CP violation we choose Am? = 10“"^eV^ and AM^ = 
10“^eV^, and take S12 = 1/2, S23 = l/\/2, and S13 = a/OT, as done by Arafune et ah 
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1^. Then, the RHS of (^) is estimated to be 0.39x10“^. But, we should note that the 
parameters we picked up are not those which maximize AP 3 . 

The above discussions on relative and absolute magnitudes of CP violations can be 
conhrmed by the precise computation of AP3 as we did for AP2. We plot in Fig. 3 AP3(z/^ — 
Ue) and > i>T-) as a function of S 13 . We see that the matter effect contamination is 

small in both channels and AP3 ~ 0.5 for L = 250 km. In Fig. 4 we plot the same as in 
Fig. 3 but for L = 730 km. In this case, however, the magnitude of matter effect can be 
comparable to the genuine CP effect, which also implies that our analytic formula (|T^ is 
not accurate. 

What are the required statistics for the experiments designed for measurements of CP 
violation at 1 % level? Suppose that we design a long-baseline experiment which will produce 
N muon events at each detector in the absence of oscillation. We assume that number of 
appearance events in the dominant channel is ~ A/2. Then, the number of events in the 
minor channel may be of the order of ~ 10“^A. We have to make subtraction between 2 or 
3 detectors. If we require the accuracy of less than 10 % uncertainty in the number of events 
in the minor channel it means that N ~ 10^. Having such statistics with very narrow-band 
beam is certainly not easy to achieve, but may be possible in future Japan Hadron Project 
and the MINOS experiments. 

The potential problem, though technical, with the multiple detector difference method is 
that we have to dig a tunnel down from the earth surface to build an intermediate detector. 
We have to dig down to 1.1 km and 9.3 km for the KEK-PS—^-Superkamiokande and the 
MINOS experiments, respectively. We hope that the technical problem can be overcome at 
the stage when the neutrino experiments for measuring CP violation is on the time table. 

To establish the multiple detector difference method much more careful studies are re¬ 
quired on various aspects including beam design. Among them one of the most important 
is the effects of averaging over hnite energy width of the beam; our preliminary investiga¬ 
tion indicates the sensitivity of AP against the variation of neutrino energy. Keeping these 
problems in mind, we want to emphasize that this method may be the only practical way 
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of measuring CP violation to the accuracy of a few % level. 
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Fig. 1: We plot in (i) and (ii) AP 2 {i'^ z/g) and in (iii) and (iv) AP 2 {i'^ z^r) as a function 
of Si 2 , for the cases where only the genuine CP effect (open circles), only the matter effect 
(dotted lines), and both effects (solid lines) exist. We hxed the other mixing parameters as 
^23 = 3.0 X 10“^, S33 = 2.0 X 10“^ for the left two panels (i) and (iii) and S23 = 2.0 x 10“^, S33 
= 0.98 for the right two panels (ii) and (iv). The other parameters are hxed to be the same 
for all the case (i-iv), i.e., E=1.35 GeV, AM^ = 5 eV^, Am^ = 10“^ eV^ and 6 = 71/2 and 
L = 250 km. We also plot the approximated values for the case where only the matter (open 
squares) and the matter + CP (asterisks) effects exist except for (ii) where no appreciable 
difference between the exact and the approximated values can be seen. 
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Fig. 2 : The same as in Fig. 1 but for L = 730 km and E = 3.94 GeV. 
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Fig. 3: We plot in (i) Vg) and in (ii) Ur) as a function of 5 ^ 3 , for the 

cases where only the genuine CP effect (open circles), only the matter effect (dotted lines), 
and both effects (solid lines) exist. We hxed the other parameters as 3^2 = 0.25, S23 = 0-5, 
E=1.35 GeV, AM2 = eV^, Am^ = 10“^ eV^, 5 = tt/2 and L = 250 km. 
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